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ABSTRACT 

Aims. To explore the impact of surface inhomogeneities on stellar spectra, granulation models need to be computed. Ideally, the most 
fundamental characteristics of these models should be carefully tested before applying them to the study of more practical matters, 
such as the derivation of photospheric abundances. Our goal is to analyze the particular case of a K-dwarf. 

Methods. We construct a three-dimensional radiative-hydrodynamic model atmosphere of parameters T c (f = 4820 K, log g = 4.5, 
and solar chemical composition. Using this model and 3D spectrum synthesis, we computed a number of Fe i and Fe n line profiles. 
The observations presented in the first paper of this series were used to test the model predictions. The effects of stellar rotation and 
instrumental imperfections are carefully taken into account in the synthesis of spectral lines. 

Results. The theoretical line profiles show the typical signatures of granulation: the lines are asymmetric, with their bisectors having 
a characteristic C-shape and their core wavelengths shifted with respect to their laboratory values. The line bisectors span from about 
10 to 250 m s _1 , depending on line strength, with the stronger features showing larger span. The corresponding core wavelength 
shifts range from about -200 m s - ' for the weak Fe i lines to almost +100 m s in the strong Fe i features. Based on observational 
results for the Sun, we argue that there should be no core wavelength shift for Fe i lines of EW > 100 mA. The cores of the strongest 
lines show contributions from the uncertain top layers of the model, where non-LTE effects and the presence of the chromosphere, 
which are important in real stars, are not accounted for. The Fe n lines suffer from stronger granulation effects due to their deeper 
formation depth which makes them experience stronger temperature and velocity contrasts. For example, the core wavelength shifts of 
the weakest Fe n lines are about -600 ms" 1 . The comparison of model predictions to observed Fe i line bisectors and core wavelength 
shifts for our reference star, HIP 86400, shows excellent agreement, with the exception of the core wavelength shifts of the strongest 
features, for which we suspect inaccurate theoretical values. Since this limitation does not affect the predicted line equivalent widths 
significantly, we consider our 3D model validated for photospheric abundance work. 

Key words, stars: atmospheres - stars: late-type - sun: granulation 

1. Introduction puter power is allowing us to solve this problem numeri- 
cally using only the basic laws of hydrodynamics, includ- 

Solar gra nulation is the visible manifes tation of surface convec- i ng radiative transfer in the energy equation. A better un- 

tion (e.g.,|Braxet^lHMlMulldlJJ99j). When rising gas from derstanding of solar and stellar g ranulation is thus possible 

the convective zone reaches the photosphere, it cools down by through numerical simulations (e.g..lNordlund & Dravinslll99' 

radiation losses. Since the opacity decreases steeply with de- i Frevtag et all Il996t IStein & Nordlundl 119981: lAsplund et 

creasing gas temperature, the temperature of the photosphere 19991: iLudwig et alfll999t iRobinson et al 1 120031: 1 VogleJ 12004 

drops abruptly as hydrogen recombines and the gas becomes Asplundj 20051 ICollet et al.1 120061 120071 iTrampedachl 12007 



suddenly transparent. The upfiowing material in the granules de- Nordlund et al. 2008). 
celerates as it enters the convectively stable region due to the 

higher pressures above and flows sideways until encountering Theoretical models can be computed at will, but using them 

other granules before being accelerated downwards due to it now ls meaningless until we are reassured that they include the rele- 

being overdense and having negative buoyancy. Since convective vant P h y slcs behlnd the phenomenon they intend to represent, 

envelopes are present in other cool stars (spectral type F and later The validation of theoretical models, i.e., verifying that their 

on the main sequence, G and later for giants and supergiants), it most fundamental predictions are supported by targeted obser- 

is expected that they also experience surface granulation. vatlons ' 18 necessary before exploring the impact of their effects 

on other matters, such as the inferred chemical compositions in 



the case of stellar atmospheres. 



The physics of radiatively driven convection is complex; 
it is three-dimensional, time-dependent, non-local, and non- 
linear. Parameterized models of granulation may require a Stellar disks are usually unresolved. The exceptions are 
large number of free parameters whose physical interpreta- nearby red supergiants, for which interferometric observations 
tion is obscure. Furthermore, they do not guarantee a unique with present-day technology are in principle capable of revealing 
solution. On the other hand, the continuous growth of com- surface structures related to convection (e.g., Gilli land & Dupred 
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1996: IChiavassall2008h . In any case, there are no near-future 
prospects for direct observations of stellar granulation. 

Fortunately, granulation signatures are also present in stel- 
lar spectra: net line wavelength shifts and asymmetries. Very 
careful observations are required to dete c t the effects of gran- 
ulation i n this m anner dDravin 1 Il987bt iRarmrez etal.1 2008, 
hereafter Paper I). High quality data have been used by sev- 
eral authors to show that stars with convective envelopes ex- 
perience granulation with different degr e es of intensity and ve- 
locity contrasts (e.g lGravl ll982l 120051: lDravins|[l987al 120081: 
Alle nde Prieto et al] 120021) . Other v elocity fields may become 



prominent in more evolved stars dGrav et alj 120081) . Indeed, 
iGrav & Nagell d!989l) have defined a granulation boundary on 
the HR diagram based on observations of spectral line asymme- 
tries in hot and cool stars. They find that the shape and mag- 
nitude of the line asymmetries in hot stars (A-type and hotter 
on the main sequence, F-type and hotter for the supergiants) are 
qualitatively different from those that correspond to solar gran- 
ulation and argue that they cannot be related to deep envelope 
convection. 

Spectrum synthesis using classical, hydrostatic model atmo- 
spheres predicts symmetric lines whose core wavelengths corre- 
spond exactly to the input laboratory wavelengths. In real stel- 
lar atmospheres that experience granulation, line profiles com- 
ing from bright, hot granules are blueshifted because the mate- 
rial is rising and therefore approaching the observer while line 
profiles coming from the dark, cool intergranular lanes are red- 
shifted. The continuum level of the upflow line-profiles is higher 
than that of the line-profiles associated with the cool downflows. 
Thus, if the granulation pattern is unresolved, an observed line 
profile is dominated by the blueshifted component. The profile is 
then asymmetric and its o bserved core wavelength is blueshifted 
(e.g jDravins et alT l981). The line asymmetry can be quantified 
by the line bisector, defined as the location of the midpoints of 
horizontal segments joining the blue and red wings of the lines, 
and the core wavelength shifts due to granulation are often re- 
ferred to as convective blueshifts. 

The goal of this series of papers is to understand the surface 
inhomogeneities (granulation) present in K-dwarf stars and ex- 
plore their impact on the determination of fundamental parame- 
ters and chemical compositions derived from the star's spectral 
energ y distribu tion. In t he first part of the series (Ramire z et alj 
2008, hereafter Paper I), we showed that, with very careful ob- 
serving strategies, it is possible to not only detect but also accu- 
rately quantify the effects of granulation on line profiles in the 
spectra of a few bright K-type dwarf stars. A theoretical gran- 
ulation model has been calculated and it is described in detail 
in this paper, along with the tests made to ensure that its most 
fundamental predictions, namely their impact on the shapes of 
absorption line profiles (specifically those due to Fe i), are valid. 
This validation will allow us to explore with confidence the so- 
called "3D effects" on the determination of chemical composi- 
tions and fund amental parameters of K-dw arfs in the final part 
of this series (Ramirez et al. 2009, hereafter Paper III). 

2. The 3D model 

Assuming local thermodynamic equilibrium (LTE), we com- 
puted a three-dimensional radiative-h ydrodynamic mode l atmo - 
sphere using the methods described in lStein & Nordlundl (1998). 
A plane-parallel box of the stellar envelope was modeled by 
solving the fluid dynamics equations of mass continuity, mo- 
ment um conservation, and e nergy conservation (e.g., Eqs. 1- 
3 in IStein & Nordlu nd 1998). The effect of temperature enters 



through the equation of state and by the radiative heating rate, 
which is deter mined by solving the ra diative transfer equation 
(e.g., Eq. 7 in IStein & Nordlundll998l) . 

In the cal culation of the 3D m odel, the equation of state 
(as given by iMihalas et alj Il988l) . continuum opacities, and 
source functions were obtained by interpolation in the tables 
included in the updated MA RCS stellar atmosphere package 
dGustafsson et aljl975l l2008). The ad opted line opaci ties (opac- 
ity distribution functions) are those bv lKuruczl dl993albl) . 

Surface gravity, chemical composition, and entropy of the 
gas entering the simulation box are the fundamental parameters 
of the 3D model calculation. The effective temperature is not 
an input parameter as in ID models but it is rather adjusted by 
modifying the state of the gas entering the bottom of the compu- 
tational domain until it reaches an average similar to the target 
value of 4820 K. The other input parameters are logg = 4.5 and 
[Fe/H] = 0.0. The s olar abundances adopted are those given by 
iGrevesse & Sauvall(ll998l ). 

The model consists of 150 x 150 x 82 grid points represent- 
ing a rectangular box of the following geometrical dimensions: 
4.7 x 4.7 x 3.2 Mm (1 Mm = 10 6 m). The simulation box has pe- 
riodic horizontal boundaries, gas is allowed to escape at the top, 
while the density and energy of the incoming gas at the bottom 
boundary are adjusted to conserve the entropy. After relaxation, 
the simulation was run for about one hour of stellar time, from 
which 100 snapshots separated by 40 s were extracted and em- 
ployed for the line formation calculations presented in this study 
(the actual hydrodynamical time-step was approximately 0.2 s). 

The time evolution of the emergent intensity that would be 
observed on the surface of this model star is illustrated in Fig.Q] 
(in Sect. [3] we provide details on the calculation of the emer- 
gent intensity from the 3D model). The pattern of granules and 
intergranular lanes resembles that of the Sun but on a smaller ge- 
ometrical scale. The correlation between the emergent intensity 
of a given snapshot and those of the subsequent ones decreases 
smoothly and vanishes for images separated by about 10 min- 
utes, a number comparable to that obtained for the solar simula- 
tion. Thus, the lifetime of granules in our K-dwarf simulation is 
similar to that of the Sun. 

Fig. [2] shows the temperature structure of our model in the 
upper 1 Mm layers (this is the portion that we use for the spec- 
tral line calculations, as explained in Sect. [3]!. The rapid decrease 
of temperature with atmospheric height just below the visible 
surface (depth = 0){3 attributed to the feedback between radia- 
tion losses and decreasing opacity, is the most prominent feature. 
There, across a layer of only about 50 km, the gas temperature 
drops about 3000 K. The apparent discontinuity of the tempera- 
ture structure at depth ^0.1 Mm is due to this effect. In the re- 
gions of line formation (depth < 0), the temperature fluctuations 
have extremes that differ by about 1000K. 

Also shown in Fig. [2] are one-dimensiona l models com- 
puted with the ATLAS Ikuruczl I1979L Il993al) and MARCS 
(IGustafsson et al.lll975l) codesQ It is clear that, in the regions 
of line-formation, ID and 3D models show very similar temper- 



1 The "depth" variable corresponds to the geometrical depth, increas- 
ing inward, and it is zero at the plane where the spatially-averaged 
Rosseland optical depth along the vertical direction is equal to one. 

2 The references given here correspond to those that describe 
the basic properties of the models. Several improvements have 
been made since their publication. The models presented in Fig. [2] 
are those available at: http://kurucz.harvard.edu/grids.html 
(the "odfhew" version ) and http://marcs.astro.uu.se/ (see also 
IGustafsson et alj|2008l for more details on the updated MARCS mod- 
els). 
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Fig. 1. Snapshots of the emergent intensity (continuum) at 6151 A predicted by our 3D K-dwarf model atmosphere (details on the 
spectrum synthesis are given in Sect. |3J- Each panel shows an area of 4.7 x 4.7 Mm 2 which corresponds to the dimensions of a 
horizontal cross-section of the simulation box. The time, in minutes, is given on top of each panel. Only the last 22.7 minutes of the 
simulation are shown. 
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Fig. 2. Temperature structure of our 3D K-dwarf model atmo- 
sphere near the visible surface (the full simulation extends much 
deeper than shown here). The color code represents the number 
of grid points at a given temperature and depth, with red being 
the largest and dark blue the lowest number. The dotted green 
line and dashed black line correspond to temperature structures 
of ID model atmospheres (see legend on the top left). The depth 
scale is set to zero at the geometrical plane where the mean 
Rosseland optical depth is equal to 1; i.e., at the "visible" sur- 
face. 



ature structures if, for the latter, the average value of the tem- 
perature is obtained at each depth. We must warn, however, that 



it is because of the temperature (and velocity) fluctuations that 
the 3D model is more realistic and that this "similarity" with ID 
models should not be used to argue that ID and 3D models of K- 
dwarf photospheres are equivalent. The reason why the average 
3D and ID temperature structures are similar is that 3D mod- 
els are nearly in radiativ e equilibrium at solar metallicity (e.g., 
IStein & N ordlund 1998), which is not the case for metal-poor 
stars (e.g jAsplund et al.lll999l) . 

A snapshot of the temperature structure of our 3D model is 
shown in Fig. [3] Clearly, the strongest temperature contrasts oc- 
cur well below the visible surface . This phenomenon was re- 
ferred to as "hidden" granulation by Nordlu nd & Dravinsl (l990) 
in their study of a granulation model for a Cen B (spectral type 
Kl V) but note that the intensity fluctuations are still clearly vis- 
ible on the surface of this cool star model, as illustrated by our 
Fig.ffl 

The temperature and velocity fields near the geometrical sur- 
face of zero depth are shown in Fig.|4]for a given snapshot of the 
simulation (more precisely, this corresponds to a layer located 
50 km below the visible surface). Hot granules have the largest 
upward velocities and appear to be expanding from their cen- 
ter. The horizontal velocity field converges toward intergranular 
lanes, which are cooler and have the largest downward veloci- 
ties. This correlation between the vertical velocity and tempera- 
ture fields is illustrated also in Fig. [5] Note that the correlation is 
tighter at a given optical depth than geometrical depth. Similar 
results are predicted for the density and vertical velocity fields; 
the downflows contain higher density gas. Thus, the predicted 
granulation of the 3D K-dwarf model has many characteristics 
of the well-observed solar granulation. 
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Fig. 3. Temperature structure of a snapshot of our 3D K-dwarf 
model atmosphere near the visible surface (depth = 0). The 
cross-section on the front wall shows the overall behavior of 
the granulation temperature field, which becomes cooler in the 
higher layers but is inhomogeneous at any given depth. Two hor- 
izontal cuts are shown; one at depth = -0.15 Mm (predominant 
colors are green and blue) and another one at depth = (predom- 
inant color is yellow). The granulation pattern is clearly seen at 
depth = -0.15 Mm. At the visible surface, the inhomogeneities 
persist but to a lesser extent and the granulation pattern is much 
weaker. 




12 3 4 

X (Mm) 

Fig. 4. Temperature and velocity fields of a snapshot of our 3D 
K-dwarf model atmosphere at the geometrical depth correspond- 
ing to 50 km below the visible surface. The velocity field is 
superimposed on the temperature field and is represented with 
solid lines (contours tracking the vertical component) and ar- 
rows (horizontal component). In this snapshot, the temperature 
extremes differ by about 2500 K (cf. Fig. [2]) while the vertical 
velocities range from about -4 to +4 kms 1 and the horizontal 
component from to 5.6 km s with a mean value of 2. 1 km s . 
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Fig. 5. Top panel: correlation between temperature and vertical 
velocity of the 3D K-dwarf model atmosphere at the geometrical 
depth corresponding to 50 km below the visible surface, where 
the average optical depth is < logr >= 0.22. Bottom panel: as 
in the top panel but for local optical depth log r = 0.22. 



The fluctuations of the temperature, velocity, and density 
fields in the 3D simulation as a function of depth are quantified 
by their root mean square (RMS) values. The maximum RMS 
values of all these parameters occur at a depth of about +0. 1 Mm 
(Fig.[6]l. In this layer, the RMS values are about 1200 K, 2 km s _1 
(vertical component only), and 0.8 x 10~ 7 gem -3 . In the regions 
of line formation (depth < 0), the RMS temperature is roughly 
constant at about 200 K while the RMS vertical velocity de- 
creases with height, from about 1.3 kms -1 at the visible surface 
to 0.5 km s _I near the top of the simulation box. 

Similar 3D simula tions for the Sun dStein & Nordlundl 19981: 
lAsplund et al.]l2000bl 7W = 5777 K, log g = 4.44, [Fe/H] = 0) 
and Procvon dAllende Prieto et alj|2002l T„« = 6500 K, logg = 
4.0, [Fe/H] = -0.05) are available for comparison. The tem- 
perature and velocity contrasts are stronger in the Procyon than 
in the solar model. At the visible surface (depth = 0), for ex- 
ample, the RMS temperature in Procyon is about 8% compared 
to 4.5% in the solar case. In the K-dwarf model, the corre- 
sponding value is even lower, about 4.0%. The peak RMS ve- 
locities of the Procyon, Sun, and K-dwarf model are 5, 4, and 
2 kms -1 , respectively. Thus, there is clearly an effective tem- 
perature dependence such that hotter stars experience stronger 
temp erature and velocity cont rasts. This fact has been explained 
by iNordlund & Dravinsl d 1990t> . Basically, the convective flux 
scales with the total flux so that cooler stars transport less con- 
vective energy than hotter stars. In addition, cooler stars are 
denser and therefore lower velocities are enough to transport the 
convective flux. Note also that Procyon has a lower surface grav- 
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Fig. 6. Root mean square values of the temperature, velocity, and 
density fields according to our 3D model atmosphere (all snap- 
shots included). In the middle panel, the solid line corresponds 
to the vertical component of the velocity while the dotted line 
corresponds to the horizontal component. 



ity, which makes its photosphere less dense, thus enhancing the 
effect. 

The granulation pattern is driven by the sharp decline in the 
gas temperature due to radiation losses and the subsequent low- 
ering of the extremely temperature dependent continuous opac- 
ity. In the Procyon model, this occurs very close to the visible 
surface whereas it occurs a few tens to hundreds of km below 
it in the cooler models (about 100 km in our K-dwarf case; see 
Fig.[2]i. Thus, in addition to having stronger temperature and ve- 
locity contrasts, hotter 3D models such as for Procyon have their 
granulation pattern visible o n their surfaces, a phenomen on la- 
beled "naked granulation" bv lNordlund & Dravinsl (fl990). 

Although in the cooler models, such as our K-dwarf simu- 
lation, the strongest temperature contrasts occur below the visi- 
ble surface, the velocity fields associated with them reach larger 
heights, in particular the regions of line formation. For example, 
Fig. |6]clearly shows that that RMS velocity does not vanish up 
to 400 km above the visible surface. We therefore expect to see 
the granulation effects on absorption lines in K-dwarf spectra. 



3. Spectrum synthesis in 3D 

Single spect ral fines were comput ed with the 3D synthesis code 
"lte" (e.g jAsplund et a l. 2000b), which is described below. In 
Appendix lAl we compare some of these calculations with those 
obtained using a differe nt 3D spectrum synthesis program, Asset 
dKoesterke et al.ll2008l). which w e will use in the third part of 
this series (Ramirez et al. 2009) for the synthesis of continua, 



molecular features, and wide spectral regions. As we show in 
Appendix [A] the agreement between the results obtained with 
the two codes is excellent. 

The transfer equation was solved along several rays through- 
out the upper 0.8 Mm of the simulation box. The contribution 
from deeper layers to the emergent flux is negligible, thus jus- 
tifying the use of only the upper layers in the line calculations. 
To improve the vertical sampling, which is necessary due to the 
rapid change of the local physical parameters in this region and 
the sensitivity of the line opacities to these quantities, the upper 
0.8 Mm layers of the simulation box were interpolated to a grid 
of 82 depth points. In the horizontal direction, of the original 
150 x 150 grid, a coarser 50 x 50 grid was adopted by eliminat- 
ing two of every three grid points along the x and y axes, thus 
keeping the same geometrical extent. This has the advantage of 
significantly reducing c omputing time while k eeping the results 
essentially unchanged ( Aspl und et al.ll2"0~00al) . The model used 
for the line calculations is therefore a 50 x 50 x 82 grid. 

The time interval between the snapshots of the 3D model that 
were saved from the original simulation is about 40 seconds. For 
the spectral line calculations, we used only every other saved 
snapshot, as tests showed that using all of them produced es- 
sentially the same results. The maximum difference in flux from 
this experiment was from less than 0.01 % for weak lines to about 
0.04% for the strongest lines. Tests using one every four and one 
every eight saved snapshots showed that this maximum flux dif- 
ference increased by factors of only two (i.e., about 0.08%) and 
three (0. 12%), respectively. In addition to taking into account the 
time dependence of the granulation phenomenon, the use of sev- 
eral snapshots can be interpreted as the equivalent of increased 
spatial coverage. 

We solved the equation of radiative transfer using 8 polar 
(ji = cos 6 in the standard stellar atmosphere notation) and 8 az- 
imuthal (<p) angles. The number of frequencies adopted in each 
case was 71, with a fine frequency spacing corresponding to 
0.4 kms~'. For strong lines, this frequency sampling was not 
enough to reach the full extent of the wings and we therefore 
repeated the calculation using 7 1 frequencies with a spacing of 
1.5 kms" 1 . We did this instead of calculating the whole wide 
spectral range with fine spacing to save computing time. The 
two results were then merged into single line profiles. 

Similarly to the model atmosphere computation, LTE level 
populations and ionization fractions were adopted in the line 
calculations. Look-up tables for the equation of state and opac- 
ity for the wavelength region in which the problem spectral line 
is located were calculated before solving the radiative transfer. 
Interpolation from these tables during the spectral line synthesis 
instead of on-the-fly calculations reduced dramatically the com- 
puting time. 

Each spectral line was computed for 3 different values of 
the oscillator strength (log gf and log gf ± 0.5). Note that this 
is equivalent to computing each spectral line for three different 
abundances, namely Ap e and Ap e ± 0.5 dex (Ax represents the 
number density of X relative to hydrogen in a logarithmic scale 
where the hydrogen abundance is defined as Ah = 12). The 
adopted iron abunda nce was solar, more precisely, Ap e = 7.45 
dAsplund et alj|2005l) . 

We calcul ated line profil e s for the 119 Fe i and 13 Fe n 
lines used by Ra mirez et al.l d2007l) . which are representative 
of the line-lists adopted in many FGK stellar abundance stud- 
ies. The sources of atomic data are reliable; transition probabil- 
ities have been accurately measured in the laboratory, often by 
more than one group and averaged after checking for consis- 
tency between groups, while the collisional broadening param- 



6 



Ramirez et al.: Granulation in K-dwarfs. II. 



Table 1. Strong Fe i lines added to the line list of iRarrurez et alj 
d2007l) for the 3D calculation of synthetic line profiles. 



Wavelength 


l°gg/ a 


EP 


<r b 


a b 


(A) 




(eV) 


(a.u.) 




5741.85 


-1.67 


4.26 


725 


0.232 


6252.56 


-1.77 


2.40 


326 


0.245 


6335.32 


-2.18 


2.20 


275 


0.261 


6411.65 


-0.72 


3.65 


820 


0.247 


6430.85 


-1.95 


2.18 


272 


0.257 


6677.99 


-1.42 


2.69 


313 


0.268 



Transition probabilitie s are from the laborat ory measurements of 
the Oxford group (e.g., Blackw ell et al.ll 197^) . 
The last t wo columns are the va n der Waals damping constants ob- 
tained by Barklem et aTJ d2000h : a is the broadening cross section 
for an atom-perturber relative velocity vo = 10 6 cm s _I , given in 
atomic units (1 a.u. = 2.8 x 10~ 17 cm 2 ), and a the velocity param- 
eter, which is related to the temperature dependence of the cross- 
section (a oc r ( '-" )/2 ). 



eters, in particular the van der Waals damping constants, have 
be en obtained from rec ent theoretical calculations (see Sect. 4.2 
in Rarmre zet al.ll2007l for details and references). To improve 
the line strength coverage, we added 6 very strong Fe i lines to 
this list. The atomic data adopted for these strong lines are given 
in TableQ] 

The Doppler shifts introduced by the velocity field of the 
simulation were taken into account in the synthesis of the ab- 
sorption line profiles, as well as the thermal and collisional 
broadening. The microturbulence and macroturbulence parame- 
ters are not needed in 3D spectrum synthesis and were therefore 
not used. 

4. Disk-averaged line profiles 

4.1. Calculation of flux profiles 

The calculation of flux from intensities is straightforward, and 
very accurate provided enough rays (/j. and <p angles) are included 
in the integrations. The rotationally broadened intensity is 



/(Av,/i, (p, V sin z) 



_L f 2 " 

2?r Jo 



/(Av - V sin/sin6>cos0',yu, <p)d<p' 

(1) 



where V sin / is the projected rotational velocity of the star, and 
(/>' the latitude angle on the stellar disk. The frequency depen- 
dence of the intensity has been replaced here with a velocity 
variable Av = cAv/v. The disk averaged line profile is then 



F(Av, Vsinz) 



Jr-2n p\ 
Jo 



/( Av, fi, <p, V sin z) fidfidip . 



(2) 



4.2. Bisectors and wavelength shifts of V sin i = line profiles 

We used Eqs. [TJ and [2] to determine disk-averaged line profiles 
adopting V sin i = 0. The bisectors of these line profiles (Fe i 
lines only) are shown in Fig. [7] where they have been grouped 
according to line-depth and excitation potential (EP). 

According to the 3D model, the line bisector spat£| ranges 
from about 250 m s _1 for the strongest Fe i lines to about 



1 .0 
0.8 
0.6 
0.4 
0.2 

1 .0 
0.8 
0.6 



en u - 



0.2 

1 .0 

0.8 
0.6 
0.4 
0.2 

1.0 
0.8 

0.6 
0.4 
0.2 
















m 










EP<3.5 


ev : 






EP>3.5 


eV " 



-200 



-100 
Velocity (m 



100 



200 



3 We define the "span" as the difference in velocity between the red- 
dest and bluest points of the line bisector. 



Fig. 7. Bisectors of the disk-averaged Fe i line profiles calcu- 
lated with our 3D model. Four groups of lines, sorted according 
to their line-depth, are shown in the different panels. The line 
bisectors have also been divided into two groups according to 
the value of their excitation potential (see legend on the bottom 
panel). 



10 m s for the weakest features. The line bisectors show the 
characteristic C-shape attributed to granulation. The core wave- 
lengths of the weak Fe i features are shifted by up to about 
-200 m s _1 . The blueshift decreases for stronger lines and it be- 
comes nearly zero for features of residual core flux around 0.3. 
Interestingly, the 3D model predicts a core wavelength redshift 
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Fig. 8. Core wavelength shifts predicted by our K-dwarf 3D 
model atmosphere. The Fe i lines have been grouped into low 
excitation potential (EP < 3.5 eV, black crosses) and high exci- 
tation potential (EP > 3.5 eV, gray diamonds) lines. 



for the strongest features, which amounts up to 100 m s _1 . This 
result is discussed in more detail below. The core wavelength 
shifts are also given as a function of the line equivalent width 
(EW) in Fig. 1 

The large dispersion present between bisectors of different 
lines that have similar line strength, as shown in Fig. [7] is mostly 
due to the EP dependence of the core wavelength shifts. Note, 
for example, that lines of EW 50 mA have a shift that is about 
30 m s _1 higher for the low EP lines compared to those that 
correspond to the high EP lines, as shown in Fig. [8] Due to this 
effect, there appear to be two branches in the line shift vs. EW 
relations shown in Fig.[8j which is, however, due to the fact that 
our line selection was such that few lines of intermediate EP 
were included; most of our Fe i lines have either EP 2.5 eV or 
EP 4.5 eV. 

When the bisectors of lines of similar strength are compared 
excluding the core wavelength shifts (therefore comparing only 
their shapes), the dispersion is reduced significantly, although 
a small EP dependence for the detailed shapes of the bisectors 
remains. Qualitatively, while the shape of the lower half of the 
bisector is nearly independent of EP, for a given line strength, the 
upper half extends more towards the blue for the high EP lines. 
Therefore, for a given line depth, the bisector span is slightly 
larger for higher EP lines. 

Fig.[8]also shows the characteristic signature of granulation. 
The weakest lines, those that are formed in deep photospheric 
layers that experience the largest granulation effects, have the 
highest (in absolute value) core wavelength shifts. As the lines 
get stronger, and therefore the line formation depth decreases 
reaching higher photospheric layers with weaker granulation 
contrasts, the core wavelength shifts become lower. 

Since the correlation between velocity and temperature fields 
decreases as one reaches higher atmospheric layers (see the cyan 
dots in Fig. [9), one would expect the core wavelength shift of 
the strongest lines to converge towards zero. Even if the prop- 
erties of some high photospheric layer are such that most of the 
light emitted there comes from downfiows (as also suggested 
by Fig. |9j see green dots and blue line), we would expect the 
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Fig. 9. Velocity and temperature-velocity correlation indices as a 
function of depth in our 3D K-dwarf model atmosphere (all 100 
snapshots are shown here). The velocity index (green dots) is 
defined as ~/cr(v z ), where vT is the average value of the vertical 
velocity field at a given depth and cr(v z ) its standard deviation. 
The temperature-velocity correlation index (cyan dots) is given 
by 0.5 rC 1 £(v z - v7)(T - T), where n g is the number of model 
grid points and T the velocity field at a given depth with a mean 
value of T. The blue and red lines correspond to these indices 
averaged over all snapshots. 



redshift to decrease shortly after, provided we continue probing 
higher layers with stronger absorption lines. A hint of this effect 
is seen in Fig. [8] where at EW - 180 mA the line shift, which 
is already a redshift, seems to start decreasing for higher EW 
values. 

In their s t udy of line formation in solar granulation, 
Asplundetal. (2000b) encountered a problem for these pre- 
dicted convective redshifts (see Figs. 11 and 12 in their paper). 
For the strongest Fe i features, they are up to 200 m s higher 
than the observed ones. Since the wavelengths of the solar spec- 
trum have been accurately calibrated in an absolute sense and, in 
particular, no correlation between the errors in the wavelength 
scale and the line strength are expected in the solar atlases, this 
may suggest that the predicted core wavelength redshifts for 
the strongest lines are a numerical artifact, possi bly due to de- 
ficienc ies in the modeling of the outer boundary ( Asplund et al. 
2000a b). A test was performed to support this claim and is de- 
scribed in the next paragraph. Note, however, that non-LTE ef- 
fects and the presence of the chromosphere in real stars may 
be the dominant factors in determining the absolute core wave- 
length shifts of those lines, which would make these particular 
3D model predictions merely incomplete and not necessarily in- 
correct. In fact, observations of the solar photosphere show that a 
reversal of the granulation pattern does occur a few hundred kilo - 
meters above the visible surface (e.g., iJanssen & Cauzzii r2006). 
which would explain line core wavelength redshifts. In our K- 
dwarf model, a weak temperature-velocity correlation (of op- 
posite sign compared to that of the "normal" granulation pat- 
tern seen in continuum formation layers) is predicted at about 
-0.1 Mm (see red line in Fig . [9j> . 

A subset of Fe i features from our linelist were recomputed 
after excluding a number of upper layers (1,2, and 4) from the 
simulation. Their core wavelength shifts were then calculated 
and compared to those obtained from the original model. We 
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find that for Fe i lines of EW > 100 mA, the upper layers play 
an important role in determining the core wavelength shift. For 
example, for an EW - 130 mA line, removing 1, 2, and 4 layers 
increased those shifts by about 6, 13, and 24 m s" 1 , respectively. 
For Fe i lines of EW < 100 mA, however, the upper layers are 
unimportant regarding the core wavelength shift. The differences 
in the core wavelengths obtained with the full model and those 
derived excluding upper layers have a mean value of essentially 
zero, with an RMS scatter of about lms 4 . These lines are there- 
fore free from systematic errors due to the limited height of the 
simulation box and reliable to test the line formation of spectral 
features commonly used in stellar abundance determinations. 

For EW > 100 mA, the core wavelength redshift diminishes 
as more layers are added to the line profile calculation. It is 
tempting to conclude from this that the redshifts of the strongest 
lines predicted by the 3D model will become lower if more lay- 
ers are added to the simulation (although admittedly from this 
test alone it is clear that the correction would be insufficient to 
explain the +100 m s _1 predicted by the original simulation for 
some of the strongest Fe i features). However, the impact of the 
chromosphere on these upper layers in real stars cannot be ig- 
nored, which would make our no-chromosphere model for these 
uppermost layers unrealistic. 

Even though these predicted redshifts appear to reveal a lim- 
itation of our current 3D models, note that the line strengths 
are not affected in the same manner. For example, for the test 
described above, the corresponding change in equivalent width 
was between and 1 .5 mA. For the strongest features, removing 
4 of the upper layers resulted in a line equivalent width lower 
by only about 1%, a number that is comparable to, if not often 
smaller than, the typical errors in the measurement of line equiv- 
alent widths. Thus, the 3D model predictions for the strength of 
Fe i spectral features of EW > 100 mA are still reliable, given 
that most of the absorption occurs in layers far from the upper 
boundary of our 3D model. 

Fe n lines suffer from stronger granulation effects. Their bi- 
sectors span from about 100 to slightly more than 200 m s 
while the core wavelength blueshifts of the weak lines reach up 
to about -600 m s _1 (Fig. [TOt. For EW values between 15 and 
40 mA, which is the only EW range where our Fe i and Fe n lines 
overlap, the convective blueshifts of Fe n lines are, on average, 
about 100 m s _1 bluer (i.e., more negative) than those of Fe i 
lines. This is due, most likely, to the deeper formation depths 
of Fe ii lines. The Fe n number density increases with depth 
as higher temperatures are required to ionize the neutral iron 
atom. At these deeper layers, the intensity and velocity fields are 
strongly correlated. Although it would be ideal to look for the 
granulation signatures using Fe n lines, in real K-dwarf spectra 
only a small number of them are available for accurate line pro- 
file measurements. 

4.3. Rotationally broadened profiles 

To explore the effect of the projected rotational velocity on line 
profiles affected by granulation, we computed a few spectral 
lines using Eqs.Q] and [2] with several non-zero V sin i values and 
compared the resultant line bisectors. One typical example of 
this exercise is shown in Fig.fTTI 

Between V sin / = and 1 kms _I , the shape of the line bi- 
sector is only slightly affected by stellar rotation and it therefore 
retains most of the granulation signatures. On the other hand, 
for projected rotational velocities greater than 1 kms" 1 , the ef- 
fects of stellar rotation must be taken into account to properly 
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Fig. 11. Bisectors of the theoretical 5228.4 Fe i line computed 
for several values of the projected rotational velocity V sin i (see 
legend). 



characterize them. In general, the effect of a non-zero V sin i is 
to make the line more symmetric. As one might naively expect, 
additional symmetric broadening dilutes the line asymmetries. 

The effect of the rotational velocity on the convective 
blueshifts is also illustrated in Fig.[TT] Interestingly, although the 
line profiles become more symmetric for higher V sin i values, 
the absolute value of the convective blueshift incr eases. This is a 
confirm atio n of the "rota tion effect" discussed bv lGrav & Toner! 
d 19851) and iGravl ((1986) on the basis of their numerical exper- 
iments using two component granulation models. The distribu- 
tion of Doppler shifts for the granules in a V sin i = star ex- 
tends from a maximum blueshift value corresponding to gran- 
ules in the disk center to zero for the limb. Since the regions 
observed near the limb correspond to larger areas due to pro- 
jection effects, the distribution shows an increase towards lower 
blueshifts, peaks at a value close to the line-center convective 
blueshift, and vanishes at zero velocity. A large V sin i star redis- 
tributes the observed Doppler shifts from the granules in areas 
near the limb symmetrically about zero. Since the distribution 
does not change for areas near the disk center, the entire disk 
distribution will now have a red tail but a peak extending fur- 
ther towards the blue, thus producing the rotation effect. Note 
that this increase in the convective blueshift is an external effect 
that is not associated with the strength of the granulation inho- 
mogeneities. Rotation is not taken into account in the compu- 
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tation of the model atmosphere, only in the calculation of disk- 
averaged line profiles from the emergent intensities. 

4.4. Impact of instrumental imperfections 

The observed spectral line profiles have an additional external 
broadening due to the finite va lue of th e spectral resolution of the 
spectrograph. In our case (see lPaperll for details) we found that 
the spectral resolution (R = A/AA) was not constant but varied 
between 160,000 and 210,000 among our spectra. In addition, 
we found that the shape of the instrumental profile was slightly 
asymmetric. Here we explore the impact that these instrumental 
imperfections have on the theoretical line profiles. 

4.4.1. Variable resolution 

In order to explore the effects of a variable spectral resolution 
on the theoretical line profiles, we performed tests convolv- 
ing the disk-averaged, rotationally broadened line profiles with 
Gaussian instrumental profiles of FWHM = A/R. Although the 
actual instrumental profile of our observations was not Gaussian, 
its overall shape remained nearly constant (its FWHM, however, 
did change with time). We performed this test with Gaussians 
to avoid the complications introduced by non-symmetric pro- 
files. The latter is discussed below. We performed the test for 
two values of the projected rotational velocity: Vsin/ = and 
2.0 kms _I . The result is shown in Fig. [12] 

As it is clear from inspection of Fig.[T2l the effects of a vari- 
able resolving power (between 160,000 and 210,000) on mea- 
surements of line asymmetries are relatively small, specially if 
the profiles have already been rotationally broadened using typi- 
cal V sin i values for K-dwarfs. Lower resolution makes the lines 
more symmetric but, in our range of resolving power and typi- 
cal values of projected rotational velocity, the asymmetry is re- 
duced, on average, by only a few meters per second. 

Note that lowering the spectral resolution using Gaussian 
instrumental profiles results in an increased core wavelength 
blueshift. This is due to the fact that the lines are asymmetric 
and the asymmetry is such that the blue wing extends more than 
what is expected for a symmetric profile. In the "blurring" of the 



spectral line, this extended wing contributes more to the shallow- 
ing of the line core than the red wing, resulting in an additional 
core wavelength blueshift. 

4.4.2. Asymmetric instrumental profile 

In lPaper fl , we showed that the instrumental profile of our obser- 
vations was very similar to a Gaussian profile but not a perfect 
match. The largest departure from Gaussianity was observed on 
the blue wing, where a small bump was present and affects all 
our data. 

Fig. [13] shows how an asymmetric instrumental profile af- 
fects the shape of a line bisector as well as its core wavelength 
shift. This test was made generating a Gaussian profile and intro- 
ducing small perturbations instead of using the real instrumental 
profile for practical reasons; the shape of the actual instrumen- 
tal profile is not as smooth as a Gaussian function due to the 
n oise of th e ThAr exposures used to determine it (see Sect. 3.5 
in Paper I for details) and it is not as trivial as the Gaussian case 
to deconvolve a non-smooth profile to take into account thermal 
broadening, which represents about 5 % of the total broadening 
of the Th lines (we do use the real instrumental profile later for 
the comparison with our observations). 

Compared to the Gaussian case, our slightly asymmetric in- 
strumental profile affects both the line bisector span and core 
wavelength shift of the line profiles formed in stellar granula- 
tion. A blue bump on the profile reduces the span and shifts the 
entire line towards the red by an amount that depends on the size 
of the bump, as seen by the sequence (b),(c),(d) in Fig. [13] The 
size of each of these bumps is such that the area under the instru- 
mental profile compared to the Gaussian case is 1.6 (a), 2.8 (b), 
and 5 % (c) higher while the shifts are about 20, 40, and 70 m s _1 , 
respectively. Similarly, a red bump that produces a 2.8 % depar- 
ture from Gaussianity increases the span and introduces an extra 
blueshift of about 30 m s -1 , as shown by case (e). Note, however, 
that the instrumental profile of our observations does not show a 
red bump alone. Only for one of our observing runs, both a blue 
and a red bump were present. Panel (f) in Fig. [T3l approximately 
corresponds to that case (here the departure from Gaussianity 
corresponds to 7.8 % of covered area difference). The instrumen- 
tal profile of our observations varied with time, but the relative 
departure from Gaussianity was nearly time-independent. Cases 
(b), (c), and (f) in Fig. [T3] fairly represent o ur data, w ith (c) cor- 
responding to the bulk of them (cf. Fig. 4 in Paperll). 

Thus, we find that after convolving the predicted line profiles 
with our asymmetric instrumental profile (represented by cases 
b, c, and f in Fig. [T3l. the span of the line bisector is somewhat 
affected (it is reduced by about 10 m s -1 ) while the core wave- 
length shift changes significantly (about 20-40 m s towards 
the red) compared to the Gaussian approximation. 

4.5. Reference mean theoretical line bisectors and 
wavelength shift vs. EW relation 

For the comparison with observed data, the line profiles pre- 
dicted by our 3D model should be convolved with appropriate 
values of rotational and instrumental profiles. Given that these 
are relatively fine effects for the sample at hand and the fact that 
only one star in the sample (HIP 86400) has fundamental param- 
eters very close to those of the 3D model, we define a reference 
case of Vsin/ = 1.5 kms 1 (roughly the projected rotational 
velocity of HIP 86400, as shown in Appendix [0 and an instru- 
mental profile that corresponds to the observing run when most 
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of the data for HIP 86400 were acquired (see lPaperlh . Hereafter, 
we will refer to the theoretical line bisectors and wavelength 
shifts as those that correspond to this particular choice of ex- 
ternal broadening parameters. 

The line bisector is an excellent probe of line formation in 
stellar granulation because at each flux depth its value is sensi- 
tive to the inhomogeneities of a certain range of photospheric 
layers. We thus expect the shapes of bisectors of lines of a 
given strength to be similar. There will be, of course, differ- 
ences due to dissimilar transition properties such as oscillator 
strength, excitation potential (cf. Fig. |7J, collisional broaden- 
ing, etc. Nonetheless, grouping lines of similar strength to ob- 
tain mean the oretica l line bisectors is necessary in the context of 
this work. In iPaper il we showed that measuring individual line 
bisectors is a very difficult task and the results very noisy, but av- 
eraging bisectors of many lines of similar line strength resulted 
in very robust mean line bisectors. We thus follow the same pro- 
cedure with the theoretical results. 

The mean theoretical line bisectors for Fe i lines with resid- 
ual core flux around 0.30, 0.45, and 0.65 are shown in Fig. [14] 
The error bars shown there (shaded areas) correspond to the 1- 
<t standard deviation from the mean, which is the intrinsic dis- 
persion that we would expect from comparison with observed 
lines. For clarity, Fig.[l4]shows the line bisectors measured with 
respect to the line center. The core wavelength shifts are shown 
separately in Fig. [15] We perform this separation because the ob- 
served data is treated in that way to avoid introducing the error 
of the rest core wavelengths in the determination of the shapes 
of the mean line bisectors. 

For comparison, Fig. [14] also shows the mean line bisectors 
that are obtained from the V sin i = profiles without instrumen- 
tal broadening (dotted lines) and also those obtained from the 
Vsini = 1.5 kms _I profiles convolved with a Gaussian instru- 
mental profiles of FWHM equal to that of the real instrumental 
profile (R m 180,000; dashed lines). Clearly, the effect of the 
projected rotational velocity on the mean bisector shape is more 
important than that of the asymmetry of the instrumental profile, 
although the impact of the latter is not negligible, in particular 
for the regions near the continuum. 

The core wavelength shift vs. EW relation that results from 
the line profiles calculated with our 3D model, after convolution 
with the reference rotational and instrumental profiles, is shown 
in Fig. [15] This figure should be compared to Fig. [8] which 
shows the corresponding result for V sin i = and no instru- 
mental broadening. In addition to vary slower, the reference line 
shift vs. EW relation (Fig. [TBI is shifted upwards due to the ef- 
fects of the asymmetric instrumental profile (Sect. 14.4.21 1. which 
in our case is more important than the rotation effect that shifts 
the lines towards the blue instead (Sect. |4~3l l. 

In Fig. Q3] the lines have been separated into high (EP > 
3.5 eV) and low (EP < 3.5 eV) excitation potential lines, as was 
done in Sect. 14.21 Two independent cubic polynomial fits were 
made, one for each group, and are shown with the dotted lines 
in Fig. [15] These fits are used below for the comparison with the 
observed data. 

5. Comparison to observations 

In lPaper j, the very high spectral resolution, high signal-to-noise 
ratio spectra that we obtained for nine bright K-type dwarf stars 
were described in detail. Here those data are used to compare 
the predictions of our hydrodynamic model atmosphere and 3D 
spectrum synthesis with the observed shapes of Fe i absorp- 
tion line profiles and thus verify the accuracy of the model. In 



Figs.[16]and[T2]we show the comparison of bisectors and wave- 
length shifts, respectively, for four of our sample stars. 

5.1. Line bisectors 

In general, the basic observed shapes of the line bisectors are 
well reproduced by the model. Quantitatively, however, there are 
measurable differences. In Fig. [16] three average spectral line 
bisectors are shown in each panel. The one representing strong 
(weak) lines was obtained with features of line depth between 
0.25 and 0.40 (0.60 and 0.80), which roughly corresponds to 
EW = 80 - 140 mA (20 - 40 mA). The intermediate strength 
bisector corresponds to all other lines in between (i.e., line depth 
between 0.4 and 0.6 and EW = 40 - 80 mA). 

The mean line bisectors of HIP 96100 (r eff 5220 K) have 
clearly larger span than those predicted by the 3D model (see the 
top panel in Fig. [TBI . There are three reasons that can explain 
this. First, the effective temperature of this star is about 400 K 
hotter than the 3D model, second, its metallicity is lower than so- 
lar ([Fe/H] = -0.22), and, finally, its projected rotational veloc- 
ity is the lowest among our sample stars (V sin i 0.8 km s _1 )Q 
Note that the granulation effects, and in fact the absolute mag- 
nitude of the inhomogeneities, are argued to increase with effec- 
tive temperature (Sect. [2]i and also with decreasing metallicity, 
although probably the la tter effect is significantly smaller (e.g., 
lAllende Prieto et al J 19991) . while low projected rotational veloc- 
ities (Vsini < 1 km s -1 ) were shown not to affect the shape of 
the line bisectors significantly (Fig.fTTTl. 

In this context, it is interesting to see the case of e Eri 
(T e ff =a 5050 K), which is more than 200 K warmer than the 
3D model and yet does not show bisectors with a span signif- 
icantly larger than the 3D model predictions, in particular for 
the strong line case. Part of the reason for this is the high V sin; 
value of the star (approximately 2 km s _1 ), which, as seen in 
Fig. QT| can reduce somewhat the span of the line bisectors. 
Although it is not shown in Fig. [16] a similar result is found 
for the stars HIP 26779, HIP 88601, HIP 64797, and HIP 37349, 
which are also hotter than the 3D model by 50 to 300 K and have 
high V sin / values (between 2 and 3 km s ). Interestingly, the 
5 stars mentio ned abov e are also the most active in our sample 
(cf. Fig. B.l in [Paper j) . which suggests that chromospheric ac- 
tivity or some other activity-related effect has an important im- 
pact on the shapes of Fe i line profiles. 

HIP 86400 (r eff 4830 K, [Fe/H] = -0.05) is our refer- 
ence star because it has parameters very close to those adopted 
in the calculation of the 3D model atmosphere. A visual inspec- 
tion of the third panel in Fig. [16] shows that the agreement be- 
tween the observed and predicted mean line bisectors is excel- 
lent. The line bisector spans for the strongest Fe i features in 
our HIP 86400 spectrum are about 70 m s . Since this value 
agrees very well with the theoretical prediction, after accurately 
taking into account the effects of finite spectral resolution, de- 
tailed instrumental profile, and projected rotational velocity, we 
can conclude that the actual bisector spans in this type of star 
(solar-metallicity early K-dwarfs) are around 140 m s _1 for the 
strongest Fe i features. The exact observed value will of course 
depend on the properties of the instrument and the value of the 
projected rotational velocity of the star. 

The coolest stars in our sample, HIP 1 14622 (r eff 4740 K) 
and HIP 2331 1 (T eS 4640 K), have both relatively low Vsini 
(< 2 km s _1 ) and are about 80 and 180K cooler than our 3D 

4 For more details and references on the stellar V sin i values reported 
in this paper, see Table A. 1 in Pa per j 
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model. Their spectra are richer in spectral lines compared to 
those of our other sample stars and this is likely the reason why 
their mean bisectors show large deviations from the theoretical 
expectation near the continuum (the bottom panel of Fig. [TBI 
shows only the results for HIP 23311, which are very similar 
to those of HIP 1 14622). In fact, the mean bisectors approach 
zero velocity there, which points to the influence of numerous 
blends randomly distributed on both sides of the spectral lines. 
Nonetheless, farther from the continuum, the 3D model predic- 
tions agree quite well with the observations of these two cooler 
stars. Therefore, it is likely that the granulation contrast and con- 
vective velocities do not change significantly in dwarf stars of 
effective temperature between 4600 and 4800 K. 

A comment should be added here about the theoretical mean 
bisector of the strongest lines. Since their core wavelengths are 
shifted by up to +100 m s _I , and we concluded this to be a limi- 
tation of our model (Sect. 14.2b . the lower part of their line bisec- 
tors should be similarly affected. The strongest features that were 
used to determine the mean line bisectors shown in Fig.fTBIhave 
EW - 140 mA, which is higher than the upper limit for the reli- 
ability of our 3D model results [EW - 100 mA). The bisectors 
of EW =s 140 mA Fe i lines are well represented in the second 
panel of Fig. [7] Note that the low EP lines are the most affected 
by the incorrect redshifts predicted by our 3D model whereas the 
high EP lines have, on average, zero shift. When plotted with re- 
spect to their core wavelengths (i.e., shifted so that the lowest 
flux pixel is at zero velocity), low and high EP lines show some- 
what similar bisectors in this group of lines. Interestingly, of all 
observed lines used to construct Fig. [TBJ about 80 % have high 
EP (> 3.5 eV). Thus, this comparison between theory and obser- 
vation is mostly unaffected by the model uncertainties because, 
for this group of lines, both high and low EP lines show similar 
relative line bisectors (therefore the theoretical average strong 
line bisector shown in Fig. [16] does not change significantly if 
computed separately for each EP group) and the observations 
are dominated by high EP lines, which, when computed with 
our 3D model, are less affected, in an absolute sense, by these 
limitations. 

5.2. Wavelength shifts 

The line core wavelength shift vs. equivalent width relations for 
four of our sample stars are shown in Fig. [T7J along with the 
model predictions. In this figure, we have forced the average ob- 
ser ved velo city for the strongest lines (EW > 100 mA) to be zero 
(cf. lPaperlh . The basic signature of granulation is evident in all 
cases and it is reasonably consistent between theory and obser- 
vation, as long as we consider lines of EW < 100 mA, given 
that the theoretical predictions for the line shifts are uncertain at 
large EW (Sect g^. 

For HIP 96100, the hottest star of the sample, a visual in- 
spection of Fig. [T7] shows that the blueshifts are underestimated 
by the model below EW ^ 100 mA. This is most likely due 
to the higher effective temperature of the star compared to that 
of the model. In general, we expect the slope in this relation to 
increase with increasing temperature. 

The observed line shift vs. equivalent width relation for e Eri 
is more shallow than that of the 3D model and perhaps even 
more shallow also than those of the cooler stars shown in Fig. [17] 
In addition, there is significant scatter, even though this is the 
brightest star of the sample and we expect random errors to be 
minimum (in addition to the larger scatter of the bin-averaged 
values shown in Fig. [T7] the error bars that represent the line- 



by-line scatter are also larger compared to those of fainter stars). 
These properties are also seen in other high V sin i stars with 
clear signatures of strong chromospheric activity in our sample 
(these stars are not shown in Fig. [17] but in this context they 
behave similarly to e Eri). Interestingly, stellar activity has been 
shown to affect the shape s of the line cores significantly (e.g., 
lGravlll984t lBorreroll2008h . which explains our results, at least 
qualitatively. 

For HIP 86400, the reference star, the agreement between 
model and observation is excellent for lines of EW < 100 mA. 
The predicted slope of the relation for weak lines is remarkably 
similar to that given by the observational data. The maximum 
observed convective blueshift is about -150m s -1 . Note that the 
theoretical relation does not flatten completely below 100 mA 
but it seems to become slightly less dependent on EW for the 
weakest lines EW < 40 mA. This is also in very good agreement 
with the observations, even though blends introduce significant 
error to our measurements for the weakest lines. 

The agreement is also good for the cooler stars HIP 1 14622 
and HIP 23311 (not shown), confirming our suspicion that the 
characteristics of granulation are quite similar between these two 
stars and HIP 86400, despite the ~ 200 K difference in T eff . 

6. Conclusions 

A three dimensional radiative-hydrodynamical simulation has 
been computed for stellar parameters = 4820 K, log g = 4.5, 
and [Fe/H] = 0, using the same prescription that has been suc- 
cessful at reproducing the granulation features in the Sun. This 
K-dwarf model also predicts a granulation pattern, i.e., a correla- 
tion between the intensity (or temperature), velocity, and density 
fields, which is, however, weaker than that found for hotter stars. 
Strong temperature contrasts at a given height occur well below 
the visible surface, contrary to the case of hotter stars where this 
happens nearer their visible surface, but the velocity field still 
shows significant fluctuations in high photospheric layers. 

The 3D model calculations for typical K-dwarfs show a span 
in the theoretical line bisectors between about 10 and 250 m s _1 , 
depending on line strength, while the maximum core wavelength 
shifts (convective blueshifts) are about -200 m s _1 for the Fe i 
lines and -600 m s _1 for the Fe n lines. 

Line broadening due to the projected rotational velocity of 
stars must be taken into account when comparing the model 
predictions with observational data. The impact of the finite re- 
solving power of our observations (presented in the first part of 
this series) was carefully analyzed and taken into account be- 
fore comparing theory to observations. We explored the effect of 
both a variable FWHM and an asymmetric instrumental profile. 
The former has only a small impact on the line bisectors and 
wavelength shifts (within the range of spectral resolution of this 
work) but the latter produce significant additional shifts to the 
core wavelengths of the lines. 

Although there is a possibility that other effects, for exam- 
ple those related to stellar activity, modify the detailed shapes of 
spectral lines profiles, the agreement between the 3D model pre- 
dictions and the observations is satisfactory. If there are effects 
other than granulation affecting the line profiles, we expect their 
impact to be smaller. 

This good agreement demonstrates that our 3D model for K- 
dwarf granulation accounts for most of the factors that determine 
the detailed shapes of Fe i line profiles, as well as their strengths, 
as quantified by, for example, the line equivalent width. This, in 
turn, means that the physics invoked in these parameter-free sim- 
ulations is fairly realistic. Our model is thus adequate to explore 
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3D effect s on chemical abundance studies, which will be done in 
iPaper III 

It should be emphasized that even though our 3D model rep- 
resents significant progress relative to the ID scenario, there is 
still room for improvement. For example, the simulation pre- 
dicts convective redshifts of about 100 ms _1 for the strongest 
Fe i lines, but implications based on solar observations and their 
comparison to our K-dwarf data indicate that they should be 
nearly zero. Non-LTE effects, the chromosphere, or numerical 
artifacts associated with our 3D model calculations could be re- 
sponsible for this discrepancy. 
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Appendix A: The advanced spectrum synthesis 3D 
tool (Assd) 

A computer code for solving the radiative transfer problem 
in 3D model atmosp heres has been recently developed by 
Koesterke et alJ d2008l) and was made available for this work. 

iKoesterke et alJ (120081) "Advanced Spectrum Synthesis 3D 
Tool" (Asset) represents an improvement over previous 3D/LTE 
codes because it is not limited to the calculation of single lines, 
simple blends, or constant background opacities. Also, it prop- 
erly includes electron and Rayleigh scattering due to atomic hy- 
drogen, and uses improved (higher-order) interpolation schemes 
within the simulation grid points for the calculation of opacities 
and intensities. Asset is also capable of working with ID models, 
which allows ID vs. 3D comparisons to be made in a consistent 
manner. Furthermore, the opacity calculation is based o n a mod- 
ified version of the thoroughl y tested ID code SYNSPEC dHubenvl 
ll988UHubenv & Lanzlll995h . 

In Paper III, we will use Asset to calculate molecular and 
atomic features as well as the continuum over wide spectral win- 
dows. For the synthesis of atomic lines in this paper, h owever, 
we used the 3D code "lte" (e.g.. lAsplund et al]|2000bl) . as de- 
scribed in Sect. [3] We computed 23 of the spectral lines used in 
this paper (17 Fe i and 6 Fe n) with Asset to perform a code-to- 
code comparison. Input data were carefully checked for consis- 
tency. A typical result of the comparison of detailed line profiles 
is shown in Fig. lA.ll while the differences in the predicted equiv- 
alent widths and core wavelength shifts are shown in Fig. |A.2| 

In general, the wings of line profiles computed with lte are 
slightly deeper, in particular the blue wing. This makes the line 
bisector span slightly larger, by less than about 10 m s _I for the 
example shown in Fig. I A. 1 1 This also leads to a slightly higher 
(more negative) convective blueshift for the lte case, as shown 
in the middle panel of Fig. IA.2I The line core blueshift is higher 
(more negative, by about 15 m s~') according to lte. 

The top panel of Fig. |A. 21 shows that the equivalent widths 
obtained from the line profiles computed with both codes agree 
within 1 mA below 90 mA and the difference shows a small 
trend for the strongest lines, with the EW values from lte being 
lower by at most 2 mA. These differences can be attributed not 
only to the slightly different detailed shapes of the profiles but 
also to the continuum level computed for each case. 

Note that the line shift vs. EW relation is only slightly af- 
fected by the choice of 3D spectrum synthesis code (bottom 
panel of Fig. [A~2l ). Only a small shift upwards of 15 m s" 1 should 
be considered if Asset is used instead of lte. Since both the line 
bisectors and line shift vs. EW relation are nearly independent 
on the synthesis code adopted, our conclusions about the good 
agreement between theory and observation presented in Sect. 
are robust. 
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Appendix B: The projected rotational velocity of 
HIP 86400 

By convolving line profiles predicted by our 3D model with 
projected rotational velocity and instrumental profiles, as de- 
scribed in Sect. [4] we determined a more accurate V sin i value 
for HIP 86400. We did not perform similar calculations for the 
other sample stars because the strength and FWHM of the lines 
are very sensitive to the stellar parameters, in particular T e ff, and 
only HIP 86400 has parameters identical (within observational 
errors) to those of our 3D model atmosphere. 

We selected 14 of the "cleanest" iron lines available in the 
HIP 86400 spectrum and computed the difference between ob- 
served and predicted profiles around ±0.15 A from the line cen- 
ter using only two free parameters, the log gf value (or, equiv- 
alently, iron abundance) and Vsim. By minimizing the differ- 
ence between observed and predicted profiles, using a ^ 2 -like 
scheme, we determined the V sin / value of HIP 86400. No addi- 
tional broadening, in particular microturbulence or macroturbu- 
lence, was necessary to accurately fit the observed line profiles. 
One of our best fits to the data is illustrated in Fig. IB. II The av- 
erage of the V sin i values obtained from all lines, weighted by 
the quality of each fit, is V sin i = 1 .57 + 0.20 km s _1 (the error 
bar corresponds to the 1-cr scatter of the line-by-line values). 



x °' 8 

3 

Z 0-6 

> 

° 0.4 

ce 

0.2 
0.0 
1 .0 

x °- 8 

n b 

> 

2 OA 



0.0 
1 .0 
0.8 



3 OA 




0.0 




0.4 
0.2 
0.0 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 




(b) 



(c) 




(d) 



(e) 




(f) 



-2 2 
Av (km s~') 



4 -250 -200 -150 -100 
Velocity (m s _1 ) 



-50 



Fig. 13. Effect of an asymmetric instrumental profile (left pan- 
els) on the shape and core wavelength shift of the line bisector 
(right panels). The 5228.4 A Fe i line was used here to illustrate 
the effect. The top panel (a) shows the case of a Gaussian instru- 
mental profile. 
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Fig. 14. Solid lines: mean line bisectors obtained from the the- 
oretical line profiles computed with our 3D model and adopt- 
ing Vsin; = 1.5 kms 1 and an (asymmetric) instrumental 
profile corresponding to our observations of the reference star 
HIP 86400. The shaded areas correspond to the 1-cr dispersion 
of the theoretical line-by-line bisectors. Each bisector has been 
forced to have zero velocity at line center. The dotted lines show 
the mean theoretical line bisectors assuming V sin i = and no 
instrumental broadening while the dashed line shows the same 
mean bisectors but adopting a Gaussian instrumental profile of 
fl^ 180,000. 
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Fig. 15. Core wavelength shifts obtained from the theoretical line 
profiles computed with our 3D model and adopting Vsin/ = 
1 .5 km s~' and an (asymmetric) instrumental profile correspond- 
ing to our observations of the reference star HIP 86400. The 
Fe i lines have been grouped into low excitation potential (EP < 
3 .5 e V, black crosses) and high excitation potential (EP > 3 .5 eV, 
gray diamonds) lines. The dotted lines are cubic fits to the pre- 
dicted line shifts of the two EP groups. 
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Fig. 16. Comparison of theoretical and observed mean Fe i line 
bisectors for four K-dwarfs. These bisectors are measured with 
respect to the line core wavelengths, arbitrary shifts of -200 and 
-400 m s _I have been applied to the weaker lines for clarity. The 
theoretical bisectors are represented by the gray shaded areas; 
their extent corresponds to the l-cr intrinsic scatter that results 
from the averaging of several Fe i line profiles. The solid li nes 
with error bars correspond to the observations (see iPaper J for 
details). 
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Fig. 17. Comparison of theoretical and observed line core wave- 
length shift vs. EW relation for four K-dwarfs. The observed 
data (filled circles with error bars) represent average values in 
bins of equivalent width and have been shifted so that the aver- 
age vel ocity of the strongest lines (EW > 100 mA) is zero (see 
IPaper j for details). The solid lines correspond to the predictions 
of the 3D model (cf. Fig [15]). The two solid lines correspond to 
the relations obtained for low and high EP lines, as explained in 
Sect. 14.51 Note that the x-axis scale is logarithmic. The dotted 
line corresponds to zero line core wavelength shift. 



16 



Ramirez et al.: Granulation in K-dwarfs. II. 



CD 





- Ite 
o asset 



Fe 



6841.2 6841.4 6841.6 
Wavelength (A) 



-100 
Velocity 





6 


.< 




c 
c 


4 


CD 


2 






01 




o 





CD 






Z 


UJ 


-4 


< 






-6 






E 


50 



Fig. A.l. Left: Spectral lines synthesized with the 3D codes lte 
(solid line) and Asset (diamonds). Their difference (lte minus 
Asset) is shown at the bottom, amplified by a factor of 10. Right: 
Bisectors of the lines shown on the left panel. 
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Fig.A.2. Difference in equivalent width (top panel), line shift 
(middle panel), and EW vs. line shift relation (bottom panel) 
between the predictions of lte and Asset. Filled (open) symbols 
correspond to Fe n (Fe i) lines. In the bottom panel, the triangles 
(circles) correspond to the Asset (lte) predictions. 
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Fig. B.l. Minimum^ 2 model fit (solid line) to the observed 
6151 A line profile (diamonds) in our HIP 86400 spectrum. 
Residuals around +0.15 A of the line center are shown at 0.2 
relative flux. 



